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THE  MOTION  OF  A PLANE-PARALLEL  HEAVY  LIQUID  IN  A CHANNEL  WITH  A 
BOTTOM  WHICH  HAS  A STEP 


G.  Abduvaliyev,  Ch.  Dzhanybekov 

The  problem  of  motion  of  heavy  liquid  in  a channel  with  a 
bottom  which  has  a step  is  examined.  Such  a problem  was  first 
studied  by  N.  E.  Kochin  [1]  in  a linear  setup  and  was  resolved  for 
the  case  of  a step  of  negligible  height. 
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Pig.  1. 


Assume  that  a heavy  incompressible  liquid  moves  at  a speed  of 
Va  (Pig.  1).  The  presence  of  the  step  at  the  bottom  results  in 
a disturbance  of  the  liquid.  Our  objective  is  to  determine  the 
form  of  free  surface  EF  in  the  region  near  the  step. 


1 


— 


Employing  the  theory  of  functions  of  a complex  variable,  we 
formulate  a system  of  integral-differential  equations,  the  solution 
of  which  will  determine  the  form  of  free  surface  EP. 

Let  us  assume  that  in  the  case  of  a weightless  liquid  we  have 
a complex  potential  of  •w-V’WV'  and  a Zhukovskiy  function  of 

For  determining  function  let  us  plot  a band  of  width 

q on  plane  W onto  the  upper  half  plane  of  plane  t At  the 

same  time  the  solid  walls  PABE  and  free  surface  EF  are  plotted 
onto  half  lines  FA,  ABE  and  segment  EF  of  the  true  axis,  respec- 
tively (Fig.  2). 
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Pig.  2. 

It  is  not  difficult  to  verify  that  this  representation  is 
given  by  the  function 


WU)- if, 


(1) 


where  q is  the  flow  of  the  liquid  per  second. 


The  function 


(2) 


comes  into  consideration 


2 


where 


(2.1) 


I 


j 


is  the  Zhukovskiy  function  for  weighable  and  weightless  liquid  of 
the  scheme  being  considered,  Vi  is  the  speed  of  the  weighable 
liquid  at  infinity  upstream.  The  limiting  values  for  determining 
functions  Cfa(t)  and  ir  are  given  in  Table  1. 
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According  to  the  limiting  values  given  in  Table  1,  based  on  the 
Schwarz  integral,  for  7r*\t>o  we  have: 


E(t)  I , TRe[F.l9)'Ji-*']d9 

/TF  J y-t 


T\ 


ReCp(t) 

9-t 


Here  the  unessential  constants  equate  to  zero.  Considering  the 
constancy  of  the  speed  magnitude  vector  on  the  free  surface  (V  = 
and  function  (2),  from  the  last  equalities  we  have: 


V) 
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RUhx** 


F«)-  FJ*WJtTP\r- 


(3) 


(4) 


It  is  possible  to  determine  that  equalities  (3)  and  (4) 
satisfy  the  boundary  conditions  of  Table  1 and  give  the  desired 
Zhukovskiy  functions  in  the  case  of  weightless  liquid  and  also  with 
EnV-tV  given  for  weighable  liquid. 

On  the  free  surface  the  function  6*  V»/V , which  is  part  of 
the  expression  under  the  integral  sign  from  formula  (4),  satisfies 
the  Bernoulli  equation.  Thus  from 


we  find 


(5) 


Here 


of  the  Proude  reciprocal,  v«»  V-,  are  the  speed  vector  magnitudes 
of  weighable  liquid  at  points  with  ordinates  y«.  and  y respectively, 
y is  the  acceleration  of  gravity,  h is  the  height  of  the  free 
surface  from  the  bottom  at  infinity  upstream.  From  equality  (4). 
proceeding  to  the  limit  with  (first  designating  the  variable 

of  integration  by  , with  and  equating  the  imaginary  parts, 

we  obtain 


where 


a*  yutitTv+u-rm'-i' 


(6) 


with  , & is  constant  H>i  ). 

From  the  first  equation  of  system  (2.1)  we  have 


Fd)  •tye*  vjjfo 


or  with  i»9  , considering  that  a-ofc*iy  and  equality 


we  obtain 


dw  *2 

djmnfrr' 
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from  which  we  have  finally 


dx  coiO 

Sw-mn'-r" 
s£Xm-  . 

dw  fU’W'i  v • 


(7) 


The  parametric  equation  of  the  free  surface  and  thereby  its  form 
are  determined  from  set  of  equations  (6)  and  (7)  with  l . 

In  order  to  solve  the  obtained  set  of  nonlinear  singular  kernel 
integral-differential  equations  (6)  and  (7)  we  employ  the  Newton- 
Kantorovich  method  £2j.  As  the  initial  approximation  of  Xm^x.(%) 
and  let  us  take  the  solution  of  the  corresponding  problem 

for  weightless  liquid  and  for  arriving  at  the  first  approximation 
let  us  proceed  as  follows. 


Assuming  the  zero  approximation  is  known,  in  the  vicinity  of 
point  let  us  replace  the  nonlinear  terms  from  expressions 

(6)  and  (7)  with  the  linear  terms  relative  to  unknowns  A.li. : 


Bnti+e (y*-y):A]l  ~ln.[i*t(y*-y>)ik]-  j^tf^ ^ * ( 9 ) 

Then  on  the  basis  of  equation  (9),  formula  (6)  is  presented  as 

ft  - A ■ F.fn-  J dll—. 


(10) 


here 


f.w-  ^ 'feht'vy-'ihL  d . 


From  the  second  equation  of  set  (7)}  in  conformance  with  formulas 
(8)  and  (9),  we  obtain 
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where 
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TTTT^TT  “ f (/ 


Let  us  introduce  the  notation 


y.fg) 


(12) 
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Then  equation  (11)  has  the  form 


Here 


fj?)  6" * CvsG*  (Q,-9t+  ■ 
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(13) 


(13,1) 


(13.2) 


(13.3) 


The  obtained  linear  singular  kernel  integral-differential  equation 
is  solved  under  the  secondary  condition 

V*d-7)~0  (1*0 

at  interval  ),  . ( is  a fairly  small  number). 

Solving  equation  (13)  taking  condition  (14)  into  account,  we 
find  Vte*  and  y/£)  is  determined  in  conformance  with  formula  (12) 


(15) 


9i(*)  is  found  from  expression  (10)  by  means  of  the  computation 
of  singular  integrals 


e.cm=a(f)+F.  <*)-  ■ 


(10.1) 


After  obtaining  and  yt(f)  the  unknown  Xt(f)  Is  computed, 
according  to  the  first  equation  from  (7),  using  the  simple  quad- 
rature 
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under  the  secondary  condition 


3Ct(i-d)=X. 


Here 
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8,8ft  are  the  abscissae  of  point  B and  the  point  which  lies  between 
A and  F on  plane  t,  respectively.  Calculation  is  made  with  B = 5 
and  8,  -t,oi . 

Let  us  look  for  the  solution  of  the  problem  in  question  (13)- 
( 1 4 ) in  terms  of 

7U*).  (19) 

where  Tm(*)  are  Chebyshev  polynomials  of  the  first  kind,  CU  are 
the  coefficients  still  unknown. 


We  find 


"4r  = 2a»iTh!  (S)®  (a), 


where  'U-.JKJ  are  the  Chebyshev  polynomials  of  the  second  kind . 

Substituting' the  relationships  (19)  and  (20)  into  equation 
(13),  considering  the  following  equality  [31 


we  obtain  the  transcendental  equation 


V 


8 


from  which  for  determining  the  coefficients  , substituting  £-ju 

(£,—  Coatf  are  the  colocation  points),  we  arrive  at  the  algebraic 
system 

C»m,m  Sr,  , f/-{)  , . 

wi#  ^ f21) 

where 

C«m  - X ff./T;  (W*b[m(  4-1 1 J 4.  ^ ^ 

& = ’ (22) 

Here  ur(?)=0,  U-tf?)* o,  T,^)m4.  The  coefficients  are  determined 

from  the  set  of  algebraic  equations  (21). 

The  Newton-Kantorovich  method  used  has  a quadratic  convergence, 
the  convergence  of  the  process  depending  essentially  on  how  well 
the  Initial  approximation  has  been  selected.  The  solution  of  a 
specific  problem  for  a weightless  liquid  was  taken  as  the  initial 
approximation,  which  is  considered  most  appropriate.  Calculations 
were  conducted  on  a "Minsk-22"  computer  for  e=0,0i,£*0.i ; £ with 

and  for  f- * with  . 

Numerical  calculation  shows  that  in  the  case  jc<  * : 2 iteration 
converges  for  t*o,i  in  the  third  approximation,  and  for  £*0,2  - 
in  the  fourth  approximation.  But  if  x^^.then  the  iteration  con- 
verges for  C*(Xi  in  the  third  approximation  and  for£«Q2  in  the 
fifth  approximation. 

In  our  calculations  we  were  guided  by  the  following: 
the  process  of  iteration  was  terminated  when  the  following  con- 
ditions had  been  met 

£*0.1.%,...  . 

BIBLIOGRAPHY 

1.  Kochin,  N.  E.  Collected  works,  v.2,  19^9. 


9 


2.  Kantorovich,  L.V.  Functional  analysis  and  applied  mathematics. 
UMN,  3,  No.  6 (28),  1948 . 


3.  Beytman,  G.,  Erdeyi,  A.  Higher  transcendental  functions,  v.2, 

M.,  1966. 


10 


DISTRIBUTION  LIST 
DISTRIBUTION  DIRECT  TO  RECIPIENT 


ORGANIZATION 

microfiche 

ORGANIZATION  MICROFICHE 

A205 

DMATC 

1 

E053 

AF/INAKA 

1 

A210 

DMA  AC 

2 

E017 

AF/RDXTR-W 

1 

B344 

DIA/RDS-3C 

9 

E403 

AFSC/INA 

1 

C043 

USAMIIA 

1 

E404 

AEDC 

1 

C509 

BALLISTIC  RES  LABS 

1 

E408 

AFWL 

1 

C510 

AIR  MOBILITY  RfcD 

1 

E410 

ADTC 

1 

LAB/PI 0 

C513 

PICATINNY  ARSENAL 

1 

FTD 

C535 

AVIATION  SYS  COMD 

1 

CCN 

1 

C591 

FSTC 

5 

ASD/FTD/  NIIS 

3 

C619 

MIA  REDSTONE 

1 

NIA/PHS 

1 

D008 

NISC 

1 

NIIS 

2 

H300 

USAICE  (USAREUR) 

1 

POOS 

DOE 

1 

P050 

cia/crb/.add/sd 

1 

NAVORDSTA  (50L)  1 
NASA/KSI  1 
APIT/LD  1 
IiLL/Code  L- 389  1 


FTD-ID (RS)T-1853-78 


